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Abstract 

The question: whether quantum coherent states can sustain deco- 
herence, heating and dissipation over time scales comparable to the 
dynamical timescales of the brain neurons, is actively discussed in 
the last years. Positive answer on this question is crucial, in partic- 
ular, for consideration of brain neurons as quantum computers. This 
discussion was mainly based on theoretical arguments. In present 
paper nonlinear statistical properties of the Ventral Tegmental Area 
(VTA) of genetically depressive limbic brain are studied in vivo on the 
Flinders Sensitive Line of rats (FSL). VTA plays a key role in gener- 
ation of pleasure and in development of psychological drug addiction. 
We found that the FSL VTA (dopaminergic) neuron signals exhibit 
multifractal properties for interspike frequencies on the scales where 
healthy VTA dopaminergic neurons exhibit bursting activity. For high 
moments the observed multifractal (generalized dimensions) spectrum 
coincides with the generalized dimensions spectrum calculated for a 
spectral measure of a quantum system (so-called kicked Harper model, 
actively used as a model of quantum chaos). This observation can be 
considered as a first experimental {in vivo) indication in the favour of 
the quantum (at least partially) nature of the brain neurons activity. 



PACS. 87.19.La Neuroscience - 87.18.Sn Neural networks - 87.17.-d Cellular 
structure and processes 



1 Introduction 



Many authors have argued that consciousness can be understood as a macro- 
quantum effect. In particular, Penrose [T| proposed that this takes place in 
microtubules, the ubiquitous hollow cylinders that among other things help 
cells maintain their shapes. It has been argued that microtubules can pro- 
cess information like a cellular automaton [2] , and Penrose suggests that they 
operate as a quantum computer (see for a review [3]). On the other hand, 
a number of other authors have conjectured that environment-induced deco- 
herence will rapidly destroy the quantum macrosuperpositions in the brain 
(see for instance and references therein). Indeed, the decoherence (as 

well as dissipation and heating) is a serious obstacle to all applications ex- 
ploiting quantum coherence. Recently, considerable effort has been devoted 
to designing strategies able to counteract the undesired effects of the cou- 
pling with an external environment. Notable examples of these strategies in 
the field of quantum information are quantum error correction codes, error 
avoiding codes [U], and "parity kicks" method [7j. It is only beginning of 
such activity. It is quite possible that their exist still unknown strategies 
(and processes [E],[I]]) which can eliminate in principle any undesired effect 
of the environment. Nature can utilize these strategies (processes) in neu- 
rons. Therefore, the simple calculations of the decoherence characteristic 
times (without taking into account such corrective possibilities) cannot be 
considered final verdict. 

Methods of statistical physics are actively used for investigation of prob- 
abilistic properties of different neurons JU|-[2n]- in the present paper we will 
study the data obtained in vivo for so-called Ventral Tegmental Area (VTA). 

The Ventral Tegmental Area (VTA) is a midbrain nucleus consisting of 
the dopaminergic cells. VTA is known as a part of the limbic brain that cor- 
responds to such high brain functions as cognition, learning, rewarding and 
emotional behavior. VTA plays a key role in the generation of pleasure and 
in development of psychological drug addiction. This area is also involved in 
control of the gonadal hormones |2Ij . 

The VTA dopaminergic cells fire in irregular manner with a mean rare 
rate between 0.5 to 10 Hz. Firing patterns of the healthy VTA dopaminer- 
gic cells include single spikes and short bursts containing 3-7 spikes. It is 
believed that the bursting manner of firing is probably pulled on by reward- 
ing stimuli accepted from the glutamate neurons originating from prefrontal 
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cortex and hippocampus. Bursting activity of the normal VTA dopaminer- 
gic cells results in dopamine release in limbic brain. Dopamine release from 
VTA dopaminergic cells accompanies rewarding behavior. VTA dopamin- 
ergic cells are autoregulated via dopamine autoreceptors expressed on their 
somas. Additionally, activity of the VTA dopaminergic cells is regulated by 
serotonin, noradrenalin and acetylcholine. 

Specific nonlinear properties of the VTA dopaminergic neurons are ac- 
tively studied in recent years (see, for instance [13], [Ej an d references therein). 
Relationship between these properties and rewarding function of VTA is one 
of the topical problems. Therefore, comparative analysis of the nonlinear 
properties of the VTA signals generated by normal brains and ones gener- 
ated by brains with genetically suppressed rewarding function of VTA seems 
to be of significant interest. For this purpose we use genetically defined rat 
model of depression (Flinders Sensitive Rat Line - FSL). Moreover, the ge- 
netically defined depression of the FSL rats (which suppresses the bursting 
activity of the dopaminergic neurons) and anaesthesia used in our experiment 
(see below) allows us extract a background signal generated by the dopamin- 
ergic neurons. Using universality of the background signal we will combine 
the signals taken from 5 single dopaminergic neurons (3 FSL rats) in order 
to obtain statistically representative data set (about 10000 spikes). 

Usually, the time intervals, T, between neighboring spikes are used in or- 
der to study statistical properties of neuron signals. In the present paper we 
will use inverse interspike intervals (interspike frequencies) uj = T _1 . While 
for probability density analysis there is no advantage in such choice, for a 
moment (multifractal) analysis the interspike frequencies turned out to be 
more informative (see below). In particular, we will show that the FSL VTA 
neuron signals exhibit multifractal properties for the interspike frequencies 
on the scales where the healthy VTA dopaminergic neurons exhibit bursting 
activity. The obtained generalized dimensions spectrum corresponds to the 
multifractal Bernoulli distribution and for the high moments this spectrum 
coincides with generalized dimensions spectrum calculated for a spectral mea- 
sure of a quantum system (so-called kicked Harper model, actively used as a 
model of quantum chaos). 
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2 Experimental methods and materials 



Male Sprague-Dawley tars were used in all experiments. Animals were anaes- 
thetized with chloral hydrate (400 mg/kg., i.p.) and mounted in stereotaxic 
apparatus. The hole was drilled 4.2 mm. anterior from the interaural line 
and 1.0 mm lateral from the medial line. Extracellular recordings were pro- 
cessed by an electrode from VTA (8.0-8.6 mm dorsal from the lambda). The 
constant level of the anaesthesia was checked by EKG and chloral hydrate 
was added as necessary. Single unit recording was carried out by amplitude 
discrimination. Each recording from the single cell included at least 2000 
spike events. After each experiment, the recording site was marked by a le- 
sion caused by 15 mA DC for 10 sec. Brains were removed and stained with 
formalin before histological examination. Frozen sections were cut at 50 mm 
intervals. Microscopic examination of the sections was carried out aiming to 
verify that the electrode tip was placed in VTA. 



The interspike frequencies: u(n) = T(n)~ 1 , have obvious meaning in the case 
of simple periodic signals with constant period T. For a varying interspike 
interval T(n) (where n is number of a spike in the spike series) meaning of 
u(n) is not so clear, especially for a random-like T{n). Therefore, it is useful 
to introduce a measure based on the interspike frequencies through a moving 
average 



This measure has a simple meaning of an average frequency for r spikes. 
Then, one can try to analyze this measure on existence of scaling properties 
using its moments: 



Usually, normalized (dimensionless) moments (u^) / (w r ) p are used for this 
purpose. In our case, however, we have \i\ = directly from the definition. 
Therefore, the scaling exponents of the dimensionless moments (if existst) 
coincide with those of the moments (2). If the scaling exists one can use the 
exponents fi p in order to calculate the generalized dimensions D p |22j : 



3 Multifractal analysis 
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Figure 1 shows composite signal, u)(n), obtained from 5 singular dopamin- 
ergic (VTA) neurons of 3 FSL rats. Because the genetically defined depres- 
sion suppresses the bursting activity of the dopaminergic neurons and be- 
cause of the anaesthesia used in our experiment we can expect a background 
signal from the FSL dopaminergic neurons. In particular, we can expect an 
universality of the individual signals (i.e. a statistic similarity even between 
the signals obtained from the dopaminergic neurons which belong to different 
FSL rats). The combined signal contains of about 10000 data points. A few 
data points with relatively short (less than 0.01 s) interspike intervals were 
excluded from the combined data set (that is about 1% of the entire number 
of the data points). 

Figure 2 shows the moments (2) calculated for the data set shown in 
figure 1, for r = 3 — 13. Log- log scales are chosen for comparison with the 
scaling equation (2). The straight lines (the best fit) are drawn to indicate 
scaling (2) in these scales. 

Figure 3 shows the generalized dimensions D p (3) (calculated using the 
exponents \i v extracted from figure 2). The scales: D p versus \np/(p — 1) 
are chosen in figure 3 for comparison with the theoretical prediction for the 
multifractal Bernoulli distribution |2"3*1: 



The straight line is drawn in figure 3 to indicate the multifractal Bernoulli 
distribution asymptotic (4) with multifractal specific heat c ~ 0.41. 



4 Discussion 

The observed scaling interval: r = 3 — 13, covers the scales where bursting 
activity of the healthy dopaminergic neurons takes place (see Introduction). 
Analogous analysis of the control (healthy) rats did not reveal the multifractal 
behavior in this interval of scales. This observation can be interpreted as 
multifractality of the background signal of the dopaminergic (VTA) neurons. 
Moreover, clear tendency to the multifractal Bernoulli distribution, seen in 
figure 3 for the high moments (the straight line), indicates a fundamental 
underlying physics (kinetics) . In order to proceed further in this direction 
we also show in figure 3 values of generalized dimensions calculated for a 
spectral measure of the kicked Harper model actively used as a model of 
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quantum chaos [24J. The kicked Harper model is obtained upon quantization 
of the following map: 



Pn+i =Pn + Ksin(x n ), x n+ i = x n - Lsin(p n+ i) 

where K and L are some parameters. Canonical quantization thus leads to 
the one-period evolution operator 

L K 
Ul k = exp[— i— cos(his)] exp[— i— cos (a;)] 
h h 

where operator v = —id/dx, and h/2iT has to be considered as an effective 
Planck constant, playing a role similar to an incommensurability parameter 
in a quasi-periodic system. Quasi-energy eigenvalues and eigenvectors are 
determined by 

Vl.k^uj = exp[-27riw]'0w 

and, if we denote by ipo the state in which the system is prepared at t = 0, 
the corresponding spectral measure will be indicated by drj^ Q (uj): its sup- 
port is contained in the unit interval. The autocorrelation function is easily 
expressed in terms of the spectral measure as follows: 

C{t) = (ip \U L>K \ipo) = [ drj^uj) exp[-27riujt} (5) 

Jo 

Information on the spectral measure is thus obtained by inversion of (5). 

Multifractal analysis requires a sequence of approximations to the asymp- 
totic spectral measure p^, where is the probability attached to the z-th 
ball (of size /) covering the support of the spectral measure. The generalized 
dimensions is given by 

ZX,= hm 1 l ° gXkip) 



p 



k- 



+oo (p - 1) log l k 



where 



Xk(p) = J2Pj(h) p , h = l/(2N k + l) 

3=1 



Phase diagram of the kicked Harper model is roughly divided into three 
regions [24], characterized, respectively , by a pure point spectrum (region 
A), a purely continuous spectrum (region B), and a mixed spectrum (region 
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C). Within region A scaling features were consistently not observed. Re- 
gion B is, on the other side, characterized by good scaling properties, and 
a converging, non-trivial spectrum of generalized dimensions were obtained 
in [25- The generalized dimensions for the parameters K = 2, L = 5 (with 
iV = 6400) are shown as crosses in figure 3. 

Reported above results can be considered as an indication of existence of 
a background dopaminergic neuron signal with universal multifractal prop- 
erties for the scales where healthy dopaminergic neurons exhibit bursting 
behavior. This signal has multifractal Bernoulli distribution as high mo- 
ments asymptotic that may be related to an underlying physics (kinetics) 
of mesoscopic systems j2Sj- in particular, the multifractal properties of a 
spectral measure of the kicked Harper model are in good quantitative cor- 
respondence with the observed multifractal properties of the dopaminergic 
neuron signals. The kicked Harper model is a prominent example of quasi- 
periodically driven quantum systems with a chaotic classical analogue, that 
also can be of significant interest for the discussion mentioned in Introduction. 

The authors are grateful to K.R. Sreenivasan for cooperation and to the 
Machanaim Center (Jerusalem) for support. 
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Figure 1: Composite signal, u)(n), obtained from 5 singular dopaminergic (VTA) neurons 
of 3 FSL rats. 
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Figure 2: Logarithm of moments (2) versus lnr for the data set shown in figure f , for 
r = 3 — 13. The straight lines (the best fit) are drawn to indicate scaling (2). 
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Figure 3: Generalized dimensions D p versus \np/(p — 1) calculated using the exponents 
fj, p extracted from figure 2 (open circles) and for the kicked Harper model [2] (crosses). 
The straight line (the best fit) is drawn in figure 3 to indicate the multifractal Bernoulli 
distribution asymptotic (4). 
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